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ABSTRACT 

We present a spectroscopic redshift of z = 1.675 for the well-known multiply lensed system 
of arcs seen in the z = 0.39 cluster C10024-I-16. In contrast to earlier work, we find that the 
lensed images are accurately reproduced by a projected mass distribution which traces the 
locations of the brightest cluster ellipticals, suggesting that the most significant minima of the 
cluster potential are not fully erased. The averaged mass profile is shallow and consistent with 
predictions of recent numerical simulations. The source redshift enables us to determine an 
enclosed cluster mass of M(<100kpc/h)=l.ll it 0.03 x lO^^/i~^M0 and a mass-to-light ratio 
of M/LB(<lOOkpc/h)=32Oh(M/LB)0, after correction for passive stellar evolution. The arc 
spectrum contains many ionized absorption lines and closely resembles that of the local Wolf- 
Rayet galaxy NGC4217. Our lens model predicts a high magnification (~20) for each image 
and identifies a new pair of multiple images at a predicted redshift of z=1.3. 



Subject headings: cosmology: gravitational lensing 
clusters: individual (C10024-M6) - 



- cosmology: observations — galaxies: 
galaxies: distances and redshifts 



1. INTRODUCTION 

The z = 0.39 rich cluster C10024+16 (Zwicky 1959) 
displays one of the finest examples of gravitational 
lensing. Four clearly related images are identified 
around the tangential critical curve in HST WF/PC- 
1 data (e.g. Small etal. 1996). A further ra- 
dially directed image of the same source was later 
found in a refurbished HST WFPC-2 image by Col- 
ley etai. (1996). These arcs have been used by Colley 
etal. (1996) to construct an image of the source and 
by Tyson etal. (1998) to examine details of the mass 
distribution. 

For many years the redshift of this lensed source 
has eluded identification despite long exposures on 
large telescopes. The blue colour and lack of opti- 
cal emission lines suggests a redshift 1 <z<2 (Mel- 
lier etal. 1991). The importance of the redshift for 
lensing studies lies primarily in measuring the cen- 
tral mass and mass-to-light ratio of the lensing clus- 
ter, and in the case of C10024-I-16 these quantities 
can be measured particularly accurately for this clus- 
ter because of the symmetric arrangement of the im- 
ages. The results will depend on the source redshift 
through its effect on the ratio of lens to source separa- 
tions, dis/ds- For a nearby lens this ratio saturates 
very quickly with increasing source redshift, but for 
a more distant lens like C10024-I-16 there is a larger 
range of dis/ds and hence a larger uncertainty in the 



mass, depending on the source redshift. 

2. OBSERVATIONS 

The HST imaging data used here is that obtained 
by Colley etai. (1996). The archival data was first re- 
duced for the purposes of selecting lensed targets for 
mulislit spectroscopy. Briefiy, the images comprise 
seven exposures totaling 8400s in the F450W band 
and six exposures totaling 6600s in the F814W band. 
Images were aligned to the nearest integer pixel and 
cosmic-ray rejected to obtain an average fiux. With 
multislits on LRIS at the 10m Keck II telescope a 
15"slit was centered on a high surface-brightness fea- 
ture of the largest, and hence most magnified arc 
(upper end of arc A, see Fig 2a) to maximize the 
detection of spectral features. The total exposure 
time was 80 minutes in 0.8 arcsec seeing, using the 
300 line grating blazed at 5000 A, providing a use- 
ful wavelength range of 4500A-9500A. The resulting 
spectrum is shown in Fig 1, revealing many absorp- 
tion lines, with a redshift of z=1.675. The spectrum, 
observed at an average resolution of 6A in the rest- 
frame matches closely that of the nearby starburst 
galaxy NGC4214 (Leitherer etai. 1996), showing a 
similar continuum slope and common absorption lines 
of Sill 1527A, CIV 1550A, Fell 1608A,2344A,2382A, 
2600A, AlII 167lA, Aim 1859A. We also detect a 
foreground z=0.18 blue dwarf galaxy (Fig 2a) which 



appeared initially to be related to the lensed arcs but 
is difficult to reproduce in the lens model (Fig 2b). 

3. MODELLING THE CLUSTER LENS 

Lens models for C10024+16 have generally im- 
proved with higher quality imaging. Using ground 
based data, Kassiola etal. (1994) and Wallington 
etal. (1995) reproduced fits to the close triplet of arcs 
(A,B,C of Figure 2) but considered arc D an unlikely 
counter image. Subsequently HST images revealed 
that A,B,C and D are morphologically similar in de- 
tail (Small etal. 1996) and that a further radially 
directed arc in the cluster center, E, is another com- 
plete image of the same source (Colley etal. 1996). 
Recently a 512 parameter fit to the resolved imaging 
data for the arc system has been presented by Tyson 
etal. (1998). This solution required the inclusion of 
a number of small dark deflecting 'mascons' around 
each of the images to offset the symmetry of a domi- 
nant central potential (see Fig 2 of Tyson etal. 1998). 

To investigate the uncertainties in mass we re- 
visit the mass model using a simple approach which 
we nonetheless find sufficient to reproduce the basic 
properties of the image configuration. We start by 
assigning profiles to the brightest cluster members 
using the form advocated by NFW (Navarro, Frenk 
and White 1995), which has a characteristic scale, r^, 
and dimensionless normalization relative to the cos- 
mological critical density, 5c = Ps/Pcrit^ and allows a 
wide range of mass concentrations. Integrating the 
mass along a column, z, where r^ = (Cr^s)^ + -2^ gives 
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resulting in a defiection angle a{9) 

the image plane at position 9 = ^r'^sldh- 

Only the brightest 8 cluster members need be in- 
cluded to generate an accurate fit (all cD galaxies, see 
Fig 2a) with Vg and be unconstrained, corresponding 
to a simple vector addition of defiection fields. In- 
clusion of many fainter members produces noise as 
the fit rapidly becomes unconstrained for the small 
number of independent constraints (5 images). The 
fit is achieved with the "downhill simplex" algorithm 
(Press et al. ) by minimizing the difference between 
the model predicted locations of the three obvious 
features (HII regions) common to the 5 main images 
of the source i.e. 






a sum over all k points of all images. 

The projected (2-D) mass density contours shown 
in fig 2 are modulated by the cluster members despite 



the projection. The fit favours an overall shallow pro- 
file centered on the central tight clump of luminous 
galaxies (fig 2a, b). The two outer ellipticals are seen 
to be responsible for the largest departure from sym- 
metry generating the triple images A,B,&: C. The so- 
lution although good is not of course unique, since a 
set of discrete profiles generates a fairly smooth po- 
tential (see Fig 2b) and for this reason we do not 
need to introduce additional unknown parameters to 
describe a diffuse component. We can convert the the 
azimuthally-averaged projected model slope of ^~o.55 
in the vicinity of lOOkpc/h for comparison with an 
NFW profile. The conversion to a real space slope is 
7 = — 1 — 2^%-, or 7 = —1.26 at the critical radius. 
This is very close to the NFW expectation for mas- 
sive clusters which have a predicted slope of ~ —1.3 at 
--^lOOkpc/h (NFW, Ghingna etal. 1998) correspond- 
ing to a combination of r^ ~ 400kpc/h and 6c ~ 8000. 

4. MASS/LIGHT RATIOS 

The rest-frame luminosities of the cluster galaxies 
are converted from data numbers to an ST-magnitude 
using header information. The instrumental colours 
of the bright ellipticals are very similar with a mean of 
y^i50W — -^8i4W^=l-62 (2.818 in the AB system or 3.334 
normalised to Vega) which in the Johnson system 
B— 1=3.7 (Holzman etal. 1996) and corresponds well 
to a passively evolved old stellar population which 
at the observed redshift z = 0.39, with Zf = 3 and 
T = .01 (Bruzual & Chariot 1995) predicts colours of 
B - 7=3.825 (3.727) for n = 0.1(0 = 1) and h=0.5 
(solar metallicity). Correction to the present requires 
removal of 0.™49-0.'"52 of passive evolution. 

The arcs define a convenient radial aperture for 
comparing mass and light with a radius of 30.5 
arcsec radius intersecting the four bright tangential 
images. Using the observed lens and source red- 
shifts, we normalise the mass distribution by the ratio 
dLs/ds = 0.61 (virtually independent of curvature). 
Integrating over our model mass distribution yields 
M{e) = 1.28 X IO^^/j-^Mq for n = 0.1, correspond- 
ing to a metric radius of ~ 106kpc/h at the lens. 
Note Q enters only via the angular-diameter distance 
of the lens, so that for $7 = 1 the mass is lower by 
9% for the same angular aperture. The model mass 
is close to the Einstein mass of a symmetric lens, 
M{e < 30.5")crit = 1.33 X IO^^Mq, {n = 0.1 and 
independent of profile) as expected given the near 
circular arrangement of the images about the cluster 
which means we have confidence in the mass to an 
accuracy of less than --^2%. 

The integrated I-band light in this aperture is 
iST^ = 16.2 or Lb{9) = 3.94 x IO^/j-^Lq {n = 
0.1) after subtracting passive evolution above and 
hence, the central a mass-light ratio is, M/Lb{9) = 



324/i(M/Lb)0 at z=0 (using Mbq = 5.48). This 
is slightly higher than other lensing clusters (Kneib 
etal. 1996, Natarajan etal. 1998), but note that 
neglecting evolution reduces M/Lb by 50%, and 
should be allowed for in accurate comparisons be- 
tween clusters. This value is short of that necessary 
for closure estimated from local redshift surveys in 
the same passband estimated requires {M/LB)crit = 
15OOtIoo/i(M/L)B0 (Efstathiou, Ellis & Peterson 
1988) corresponding to a larger volume and a later 
mean galaxy type. 

The true size of the source galaxy after subtrac- 
tion of the deflection field of image radius, is found 
to be r --^0.25", ^ 20 times fainter than the tan- 
gential images , with an unlensed luminosity Mb = 
—20.85 lb 0.7 — 51og/i and an apparent magnitude 
of / ~ 24.8, typical of what may be expected for 
a galaxy at the measured redshift (see Figure 18 of 
Bouwens, Broadhurst & Silk 1998). 

5. NEW MULTIPLE IMAGES 

New multiply lensed images may be sought with 
our lens model to check and improve upon the model. 
The unknown redshift introduces an additional free 
parameter from the distance dependence in the bend 
angle. Hence to search for new multiple images we 
need only take the deflection field a{6) defined for the 
five arc system at zi = 1.67 with {dLs/ds)zx and mul- 
tiply by a scalar, / = {dLs/ds)z2/{dLs/ds)z^, mim- 
icking the effect of a change in the source redshift. 

In practice it is difficult to securely identify new im- 
ages, mainly because galaxies are similarly blue and 
numerous so that a unique identification based on 
only 2 passbands is difficult. Furthermore arcs which 
are obviously radial or tangential lie close or straddle 
the critical curves and therefore much more magnified 
than their counterimages which consequently may be 
too faint to detect. Furthermore, for a given mass 
distribution a source must lie above some minimum 
redshift to generate multiple images, corresponding 
to z>1.0 on average for sources within the einstein 
ring of C10024-I-16. Here we claim a secure iden- 
tification of one new pair of arcs as shown on fig- 
ure 2. The relative deflection is 90% of that of the 
main arcs corresponding to a predicted redshift of 
z=1.31,z=1.34,z=1.33 for $1=0.05,17=1, f^ + A=0.3 + 
0.7. In principle then a sufficiently accurate lens 
model can produce a geometric constraint on the cos- 



mological curvature, however many more multiple im- 
ages of higher redshift sources need to be identified 
to make this practical. 

6. DISCUSSION AND CONCLUSIONS 

It is clear from the above modeling that some de- 
gree of mass substructure is required in C10024-I-16, 
contrary to the conclusions of Tyson etal. (1998). 
The large mass-to-light ratio we assign to the lo- 
cation of the luminous ellipticals is far in excess of 
isolated elliptical galaxies, meaning these galaxies 
simply trace well local minima of a general poten- 
tial. This finding is consistent with the substruc- 
ture apparent in all carefully studied lensing clusters, 
notably, A2218, A370, AC114, A2390, C10939-b47 
(Kneib etal. 1996, Small etal. 1996, Abdelsalam 
etal. 1998, Pierre etal. 1996, Natarajan etal. 1997, 
Seitz etal. 1996). It is not clear if this level of sub- 
structure is in excess of N-body predictions which, 
as Ghigna etal. (1997) point out are certainly un- 
derestimates within the central 50Kpc/h where the 
problem of "overmerging" is still significant despite 
their superior dynamic range. 

Our azimuthally averaged mass profile is shallow, in 
good agreement with recent high-resolution CDM N- 
body simulations, corresponding to an NFW profile 
with Tg ~400kpc/h which is four times greater than 
the observed einstein ring radius and leads to large 
central image magnifications. Although our model 
does not explicitly incorporate a separate diffuse com- 
ponent it can be seen in Fig 2, that the sum of the 
eight profiles forms a generally smooth mass distribu- 
tion illustrating the degeneracy of this sort of model- 
ing and hence the redundancy of a separate provision 
for diffuse matter. 

We have successfully used our lens model to find 
new multiply lensed galaxies. It is clear that with 
more color information further systems will be dis- 
tinguished in an iterative procedure where the mass 
model is successively refined with the incorporation 
of the new images. Relative distance predictions can 
be made this way for comparison with measured red- 
shifts allowing, in principle, a geometric measure of 
the cosmological curvature. 

We thank Rychard Bouwens and Ben Moore for 
useful conversations. TJB acknowledges NASA grant 
AR07522.01-96A. 
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Fig. 1. — The upper curve is fluxed spectrum ofi the upper HII region of arc C (fig 2b). Many weak absorption lines are visible 
yielding an unambiguous redshift of z=1.675. Note the similarity with the local "Wolf-Rayet" galaxy (lower curve) NGC4214 
(Leitherer 1996). The sky spectrum is also shown, in the lower panel 
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Fig. 2. — The central region of the summed B+I band images is shown on the left for comparison with the model on the right. 
Note the good agreement with the 5 main images (A-E). These may be improved further by hand, incorporating masses local to 
each image, but will not alter the conclusion that the central mass distribution traces well the location of the brightest ellipticals. 
A new multiple image pair is identified by the model (i,ii) at a predicted redshift of z ~ 1.3. The contours represent surface mass 
density from (0.6 — 1.2)Ec,it, separated by 0.1. 



